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Abstract 

A general framework is presented for the renormalization of Hamiltonians 
via a similarity transformation. Divergences in the similarity flow equations 
may be handled with dimensional regularization in this approach, and the 
resulting effective Hamiltonian is finite since states well-separated in energy 
are uncoupled. Specific schemes developed several years ago by Glazek and 
Wilson and contemporaneously by Wegner correspond to particular choices 
within this framework, and the relative merits of such choices are discussed 
from this vantage point. It is shown that a scheme for the transformation of 
Hamiltonians introduced by Dyson in the early 1950's also corresponds to a 
particular choice within the similarity renormalization framework, and it is 
argued that Dyson's scheme is preferable to the others for ease of computation. 
As an example, it is shown how a logarithmically confining potential arises 
simply at second order in light-front QCD within Dyson's scheme, a result 
found previously for other similarity renormalization schemes. Steps toward 
higher order and nonperturbative calculations are outlined. In particular, a 
set of equations analogous to Dyson-Schwinger equations is developed. 
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I. INTRODUCTION 



While particle physics has been and remains to be guided by the desire to understand 
phenomena at ever shorter distances, one may nevertheless argue that the most vital chal- 
lenge for quantum field theory has little to do with the physics at Planck length scales or 
a theory of everything (TOE). This challenge is how to solve strongly interacting quantum 
field theories, which inevitably amounts to understanding in a nonperturbative way how 
a theory designed to fit the physics of a smaller scale will determine the phenomena ob- 
served at greater length scales. The prototype is Quantum Chromodynamics (QCD), but 
the problem has existed since soon after the development of the renormalization program 
for Quantum Electrodynamics, when physicists began to study strongly interacting meson 
theories; and the problem must presumably be confronted by anyone who wants to show 
that any given TOE is indeed a TOE. 

In ordinary quantum mechanics, strongly interacting systems may be handled by a variety 
of techniques based upon a Hamiltonian formulation of the theory in question. Relativistic 
quantum field theory, on the other hand, is studied almost exclusively in the Lagrangian 
formulation nowadays. Lagrangian methods are favored because they better allow one to 
keep symmetries of the theory manifest. However, this generally entails a sacrifice in calcu- 
lational flexibility, with the result that one is limited to applications for which perturbative 
or semiclassical methods are valid (or at least hoped to be valid). This paper takes the 
view that a Hamiltonian formulation will provide a better starting point for the study of 
strongly interacting quantum field theories. It was this motivation, and the difficulties in 
its realization, that eventually led Wilson to formulate his version of the renormalization 
group. 

The renormalization group is most often used in a very restricted way in quantum field 
theory. Typically, renormalization is used as a tool only to cover up infinities, to select out 
well-behaved ( "renormalizable" ) theories, and to make finite redefinitions of couplings so 
as to reorganize perturbative expansions in a better way. Renormalization transformations 
that produce more complicated interaction terms are generally avoided. This is because 
such transformations usually break at least some symmetries and complicate perturbative 
expansions. For a theory like QCD, however, perturbation theory is not powerful enough 
to calculate low-energy phenomena, for the lowest orders of a perturbative expansion look 
nothing like what finds from low-energy experiments. In such a case, it might be worthwhile 
to consider a renormalization group transformation that brings new interactions, even if such 
a transformation hides symmetries that in the usual perturbative scheme remain manifest. 
The idea is that such a renormalization group transformation can be designed so that the 
new interactions it generates will contain the most important physics. 

In fact, this is the major motivation behind the similarity transformation of Hamilton- 
ians.[[] The object of the similarity transformation is not primarily to produce a renormalized 
Hamiltonian, but rather to produce a Hamiltonian that can be dealt with in a way such that 



^ Matrix elements of the transformed Hamiltonian H' are "similar" to the original Hamiltonian H 
if H' = PHP^ 1 . If H' is to be hermitian, we need Pt = P -1 , and so the similarity transformations 
described here will all be unitary transformations. 
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even first approximations in its analysis provide a qualitatively correct description of the 
physics. That the similarity transformation can be used to renormalize a Hamiltonian that 
contains divergences is just a particular result of this more general requirement. The practi- 
cal manifestation of this basic requirement is that the similarity transformation produces a 
Hamiltonian where only states of nearby energy contribute to phenomena at a given energy 
scale. Thus Glazek and Wilson! first introduced their similarity transformation to deal with 
divergences in light-front field theory, but also with the idea of using it to help transform the 
QCD Hamiltonian into a form for which a lowest-order approximation begins to look like a 
constituent quark modelJl Wegner,§ on the other hand, independently, but at roughly the 
same time, introduced a set of Hamiltonian flow equations to deal with many-body problems 
which have no need of the renormalization of divergences. 

We shall follow the presentation of the similarity renormalization scheme given in Ref. 
[3], generalizing it in the following. The transformations of Wegner and of Glazek and Wilson 
will be seen to follow from more specific choices within our formalism. These three authors 
are generally given credit for introducing the similarity transformation to field theory. In- 
terestingly enough, however, soon after the development of renormalization theory, Dyson 
himself in a difficult, fascinating, and seemingly universally ignored series of papersu devel- 
oped what he called the intermediate representation, which results from a transformation 
that matches all the requirements of the similarity transformations developed by Glazek, 
Wilson, and Wegner. This is of more than just historical interest!: we shall see that Dyson's 
scheme has some features that lead it to be preferred over the more recent schemes. 

The outline of the paper is as follows. This section has been an introduction. In the 
next section we provide a general interaction representation formulation of the similarity 
renormalization of Hamiltonians with respect to differences in unperturbed energies. We 
show how the schemes of Wilson and Glazek and of Wegner correspond to specific choices 
in our framework and argue that one might be able to make a more manageable specific 
choice of transformation than these. In Section III we discuss the perturbative solution of 
the similarity flow equations for these schemes and conclude that the corresponding dia- 
grammatic expansion, which corresponds to that of old-fashioned time-ordered perturbation 
theory, leaves something to be desired in terms of computational simplicity. 

This leads to a detailed discussion of Dyson's intermediate representation, which we show 
corresponds to a specific choice of perturbative scheme in the similarity renormalization 
framework. The beauty of Dyson's scheme is that it may be represented diagrammatically 
in terms of generalized Feynman graphs — generalized in that amplitudes may be off energy 
shell. The difficulty comes in demonstrating this, and we dedicate some space to summarizing 
Dyson's presentation. Once the diagrammatic expansion is developed, we conclude that 
Dyson's scheme is much simpler computationally — at least for perturbative solutions - 
than those of Wilson, Glazek and Wegner. We end this long section with an example 
calculation taken from QCD. 

In Section IV we formulate the solution of the similarity flow equations in terms of 
general n-particle "dressed" vertices. The resulting equations, which look much like Dyson- 
Schwinger equations, are of second order in these vertices; and we show how a specific 
choice of similarity transformation — a modification of Wegner's scheme — can guarantee 
that these dressed vertices are joined by Feynman propagators, which therefore permits a 
Wick rotation to Euclidean space. We show how analogous equations can also be derived in 
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Dyson's scheme and that these equations are linear in the dressed vertices. Finally, we give 
a specific example of how one may make approximations to sum a set of diagrams similar 
to those summed in the Bethe-Salpeter equation. 

We shall refer to a generic theory of fermionic and vector bosonic fields throughout 
the paper. To provide a flavor of the type of calculation involved here, we give specific 
examples for QCD in light-front coordinates and light-cone gauge, but we shall only refer 
to other work that presents the Feynman rules for these particular calculations. In Section 
III, we show that the similarity-transformed Hamiltonian contains a potential that rises 
logarithmically for large quark- ant iquark separations, a result obtained previously in other 
similarity renormalization schemes. In Section IV we discuss a ladder-type summation of 
these potentials. These examples are only meant to be suggestive. More detailed applications 
of our similarity renormalization framework will be undertaken in future papers. 

Section V concludes the paper with a discussion of the ideas presented here. This last 
section is not just a summary but rather is intended to be an integral part of the paper, for 
it provides further motivation that could not properly be given before the scheme itself was 
developed. In particular, the significance of the similarity renormalization scheme for appli- 
cations in light-front field theory and the motivation for the use of light-front coordinates in 
field theory in this light are discussed at length. 

II. SIMILARITY RENORMALIZATION SCHEME 
A. General formulation 

We begin with a regularized field theoretical Hamiltonian H, derived in the usual way 
from a renormalizable relativistic Lagrangian. We write H = H + H int , with H the usual 
free Hamiltonian. Let us assume that H describes processes which at high energies may be 
well approximated by treating the interaction H int perturbatively, but that its low-energy 
processes and its bound states are essentially nonperturbative. Clearly, the problem here is 
that H is a poor first approximation for calculating low-energy phenomena. The similarity 
renormalization scheme seeks to change the original Hamiltonian into a form that is more 
suitable for calculations over a wider range of energies. This is done by introducing an energy 
scale a via a similarity transformation S a (which will always be unitary here). Roughly 
speaking, the transformed Hamiltonian H a = S a HSl will depend on a, the similarity scale, 
as follows: processes described by H which transfer an energy much greater than o will be 
absorbed into the structure of H a , so that H a will only describe processes that correspond 
to energy transfers of order a or less. 

Thus we call H the undressed Hamiltonian and H a the dressed Hamiltonian. The un- 
dressed Hamiltonian is renormalized; that is, the undressed Hamiltonian is the sum of a bare 
regularized Hamiltonian Hb and a set of counterterms He that removes any dependence in 
physical quantities upon the regulators. The structure of these counterterms must be deter- 
mined as one solves for the dressed Hamiltonian H a . In particular, the counterterms in H 
ensure that H a will remain finite as the regulators are removed. The dressed Hamiltonian 
H a thus cannot produce infinities; and so the similarity transformation S a is finite, pro- 
vided that it avoids singularities corresponding to small energy denominators. This latter 
restriction shall be discussed presently. 
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To define such a transformation, we write the unitary operator S a in terms of the anti- 
hermitian generator of infinitesmal similarity scale transformations T a : 



S a = S exp I J da'T a , 1 , (1) 

where S puts operators in order of increasing scale. Then the similarity transformation is 
given by the differential "flow" equation 

JIT 

with the boundary condition = H. Now let the eigenstates of H be |z), with H \i) = 
Ei\i); and write the matrix element of an operator O as Oij = (i\0\j). We wish to choose 
T a in such a way that 

Hoii = Gay, (3) 

where the "similarity form factor" f(x) is a smooth function such that 

(0 /(0) = 1; (4) 
(ii) /(s) ^ as x — > oo; 
(in) f*(x) = f(-x); 

and G a is a finite hermitian operator. The similarity function / has often been chosen to be 
a step function in previous work. However, this introduces nonanalyticities and in general 
should be avoided. 

It is convenient to first write H a = H + H 1 ™ 1 and transform to the interaction picture 
Oj(t) = e lH otQ e -iH t ^ go "similarity flow" equation becomes 

dH^l = _i<BM_ p>.(t), *.(*)], (5) 

where Hi a (t) = e lH otjjint e -iH t — an( ^ go ^ e undressed interaction Hamiltonian is Hj(t) = 
lim^oo H Ia (t). After Fourier transforming, (Q) and (||) become 



H Ia (u) = ^2n5(w -Ei + E 3 )H Ial] = f^G Ia (uj) (6) 

and 

= u>T Iir (u>) - [^[T IrT (u%H IrT (u-u% (7) 
do J In 

respectively. Now we are ready to construct some suitable similarity operators T a (uo). 
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B. Wilson's and Glazek's scheme 



Wilson and Glazek first introduced similarity renormalization in the context of light- 
front Hamiltonian field theory.! Their transformation was later elaborated upon in Ref. [3] . 
Essentially, their choice of T a gives 



TjJu) = — 




H Ia (u) 



(8) 



which shall be called the WG scheme henceforth.^] Note that property of / ensures that 
Tj a (u) is finite as uj — > 0, thereby avoiding possible problems with small energy denomena- 
tors. This is a great advantage of the similarity renormalization scheme. The flow equation 
for the choice (H) is 



H Ia {u) = f(i)\HM + jdo'J 



~2k 




(9) 



The term d a tHi a /(u}') on the righthand side means that when we express this as a pure 
integral equation it will contain terms of all orders in Hj a . Presumably, this makes the 
equation rather difficult to solve. 

C. Wegner's scheme 



Wegner introduced his Hamiltonian flow equation to study problems in condensed matter 
physicsfl In the present notation, his equation is 



dH a 
da 



O"' 



[[Ho, Ha] , H a ] . 



(10) 



This is a similarity transformation with the choice T a = [H , H a ]/a 3 . Actually, Wegner 
advocates that H include not just the free Hamiltonian but also the number conserving 
part of the undressed — or even better, the dressed — interacting Hamiltonian. These 
are certainly choices which one can (and, perhaps, should) make in any similarity scheme, 
although they shall not be pursued further in this work. 

Wegner's flow equation can be seen as resulting from a specific choice of the similarity 
function / in a more general scheme. This scheme, referred to hereafter as the W scheme, 
results from the following choice of the infinitesmal transformation operator: 




(11) 



T In the referred work the argument of / was 



Ej—Ej 
Ei+Ej+a- 



It is the modification to simple energy 



differences used here that makes the use of the interaction picture convenient. 
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which gives the following similarity flow equation: 



HiM = /(f) Ufa) + J da' J — [H Ia fa'),H Icrf (uj-u;'))\. (12) 



This closed integral equation looks more promising than (P), although it is still nonlinear in 
the dressed Hamiltonian. To recover Wegner's flow equation, then, one just choses f(x) = 
e~z x , which fulfills all the requirements of a similarity function. Using (j^f, the flow 
equation (|T2|) becomes 



r, / x tt ( \ 7 da ' [ du' F a {\uj + - uj') r -I 
G ItT (u) = Hfa) - / — / j — G^^u + u ),G Ia >(^u - oj ) , (13) 



where 



With Wegner's choice of /, the form factor in ( O ) is = (|a;+t<;) 2 e w2//4<J ' 2 e~ w ' 2 / 2<J ' 2 ; and the 
divergent dependence on tu, while tamed by the overall factor /(— ) in flT2|), will nevertheless 
favor higher orders in an iterative solution of Hj a (u). 

A similarity function f(x) that falls off less sharply — for example, of the form (l+x 2 )~ m , 
which goes as x~ 2m when x — >• oo — is likely to be preferred over Wegner's choice. In 
particular, for m = 1 this latter choice gives 

„ n i i /x a ,3 (^ + 2cj')(a ,2 + cj 2 ) 

F„(iw + c/, |w - w' = "7 — * T2, (15) 

(a' 2 + (fa; - u'Y) (a' 2 + {\u + u'f) 

which falls off as u~ 3 when uo — > oo. This damps far-off-diagonal contributions to Gi a {ui) 
that come from the nonlinear term in QIBp. We shall see other reasons to prefer this choice 
in the following. 



III. PERTURB ATIVE SOLUTION AND INTERMEDIATE REPRESENTATION 

A. In the WG and W schemes 

One may look for solutions to the similarity flow equations by expanding in the undressed 
interaction Hj(u) via 

G,fa) = £ / ^ • • • / ^g^J^Hfa,) ■ .-Hfa n )Hfa -co,- ... - u n ). (16) 

n=0 J 2lX J 271 

Then the functions g^^fa) = g^(oj;oJki ■■■i 0i >t), n > 1, satisfy the recursion equations 
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( n ) / \ 



m=l 



I— 



/( <j' ) (n— m) 



d 

for the WG scheme and 



1 ~~ f{j) (m-l) / \ 
_ 5W...m-lV a; mJ 



in— mj / \ 
0<r'm+l...rA u; _ U m) 



(17) 



/C" )^'m+l...n( a; ~~ w m) 



CO 



OJ, 



(18) 



fm— 1) / \ in— m) / \ 



for the W scheme, with gj® = 1 for both schemes. 

The similarity transformed Hamiltonian is then obtained by normal ordering the field 
operators in the righthand side of (|l6|) , giving a unique representation in terms of normal- 
ordered products of field operators: 



oo 

E 

m=0 ^ i 



nf d,UJi I ( m ) 



(u) : HAui) ■ ■ ■ HAum^HAu - lo x 



Urn) ■ ■ (19) 



A given normal-ordered term with coefficient g^ 1 ' in the dressed Hamiltonian will have 
contributions from an infinite number of terms with n > m in the expansion (0), with each 
contribution consisting of an integration over the product of with one or more functions 
coming from the contractions of the field operators in ([16]). Since these contractions do not 
come from time-ordered products, one has many more contraction terms than in Feynman 
perturbation theory. Indeed, the resulting diagrammatic expansion! is that of old-fashioned 
perturbation theory. The output dressed Hamiltonian is still to be solved, of course; but 
now we have a variety of nonperturbative methods, such as variational techniques, that can 
be applied. 



B. Dyson's scheme 



If the scale a is large enough, then it is reasonable to solve the similarity flow equations 
perturbatively as above; for then only high energy processes will dress Hi a , and we have 
assumed for our theory that the interactions governing such processes are weak. However, the 
proliferation of diagrams in old-fashioned perturbation theory makes the above schemes quite 
tedious. Even second-order calculations can be rough; and in light-front field theory, where 
similarity renormalization methods have been in use for several years now, little (if anything) 
has been accomplished at higher orders. Thus it is more than just interesting to note that 
Dyson long ago developed a similarity renormalization scheme which is based upon the 
evaluation of off-energy-shell Feynman diagrams. Since Dyson's intermediate representation 
was developed over several long papers, we can do it little justice here: for more details and 
the justification of many claims one should refer to the original work.i We shall oversimplify 
the presentation, change the notation, and — to a certain extent — alter the scheme itself, 
in fitting it into the present work. 
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Dyson considered the transformation Oo(t) = S , D 1 (t)(9/(t)5 , D(t) of interaction represen- 
tation operators to another representation induced by the unitary operator 

S D (t) = e* Hot S D e- iHot =r\e- i ^ dt ' Him ) , (20) 



where T is the usual time-ordering operator (acting only on the integrated times in the 
expansion of the exponential) and 

H I {t,t')=g(t-t , )H I {t') (21) 

describes a gradual switching on of the interaction: one requires g(t — t') —>■ as t' — > — oo 
and g(0) = 1. 

To be precise, Dyson introduced the function g by rescaling each occurence of the renor- 
malized coupling e R in Hjfa), including the counterterms, via e# — > g(t—ti)eR, and included 
in Hj(t,t') terms depending on time derivatives of the g(t — t$). Our modification to ( |2~T|) 
simplifies the presentation, but we should keep in mind that it is desirable to take terms of 
order (e^) n in Hj(t,t') to be proportional to g n (t — f), so that all terms in Hp of a given 
order in cr have the same strength. The modification to our presentation needed to meet 

OO , s 

this demand is straightforward: one would have #i(t,f) = £ a n g n (t - f)Hf'(1f), where 

ra=0 

Hj^ is of order (e^)™ and the dimensionless constants. 

With Sd and g chosen as above, then, low-energy processes, corresponding to large 
times, do not modify Oi>(t); whereas high-energy processes do. Indeed, since the function 
g is dimensionless, it must introduce an energy scale, which we can write explicitly via 
g(t) — > g{crt). Dyson chose the form 

oo 

g (at) = J dTG(T)e- Tat (22) 
o 

(so g is the Laplace transform of G) and showed that if one imposes the condition g'(0) = 
- that is, J dTTG(r) = — then the (Fourier-transformed) intermediate representation 
Hamiltonian obeys 

Hd(lu) oc — as oo — > oo, (23) 

demonstrating that So{t) induces a similarity transformation as defined above. Note that as 
cr — > 0, g — > 1 (which implies / dTG(T) = 1); and so in this limit Su(t) becomes the familiar 
time-evolution operator that transforms from the Heisenberg to the interaction pictures. In 
this sense, then, for finite a, Su(t) can be considered as transforming to an intermediate 
representation, which we have labelled with a D (Dyson picture) since the label I has already 
been taken. The intermediate Hamiltonian is 

H D {t) = S D \t)[H + Hj(t)}S D (t) (24) 

oo * *1 

Y, 1 * 1 j dhgiait-h)) j dt 2 g(a(t-t 2 ))--- 

o — oo 

tn-1 



J dt n g(a(t - t n )) [H z (t n ), ...[Hi(t 2 ), [#/(*i),#o + Hj 
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and the infinitesmal similarity transformation operator associated with So{t) is 



T D (t) = S D \t)±-S D (t) (25) 
= -E^y dhd a g(<r(t-t x )) j dt 2 g (a(t - t 2 )) ■ ■ ■ 

n= l — oo — oo 

in-1 

■■■J dt n g(a(t-t n ))[H I (t n ),...[H I (t 2 ),H I (t 1 )]...]. 

— oo 

Thus Dyson's transformation is an explicit perturbative representation of a similarity trans- 
formation, and the label D will always imply a a scale dependence. 

Writing H D (t) = H + Hf(t) and integrating by parts, we may express fl2l| ) as 

oo * *° ll 

H?{t) = Y, in j dt d to g{a{t-t )) J dt l9 {a{t - t{j) J dt 2 g {a{t - t 2 )) ■ ■ ■ (26) 

o — oo — oo 

tn-1 

dt n g{a{t - t n )) [Hjitn), ...[Hjih), [Hjfa), Hj(t )]]...]. 

Note that for finite a the t-dependence of the intermediate Hamiltonian is still that of an 
interaction picture operator: Hf(t) = e lHot Hfe' lHot . Now (|26"D may be expressed in the 
alternate form 

t 

Hf{t)=-a J dt g'(a(t-t )) (27) 

— oo 

oo im /An °° °° 

x V — — / dt x g{a{t -h)) ■ ■ ■ / dt m+n g(a(t -t m+n )) 

m „ =0 ™ n\ J J 

"»!'*—" — OO — OO 

xT{H I (t n+1 )...H I (t n+m )}T{H I (t )...H I (t n )} 1 

where T orders in the sense opposite to T. This latter form allows the normal-ordered 
Hamiltonian to be built up in terms of connected Feynman graphs. Let us briefly review 
Dyson's derivation. 

First, as above, we must express the intermediate Hamiltonian in terms of a sum over 
normal-ordered products of field operators. For the chronologically ordered and antichrono- 
logically ordered time products (the latter results from the Hermitian conjugate of the for- 
mer) we can use Wick's theorem. Field contractions from these are then the usual Feynman 
propagators (and their complex conjugates), and we are left with a product of two sets of 
normal-ordered operators. In normal ordering this remaining product, whose operators are 
not related chronologically, one finds field contractions that do not correspond to Feynman 
propagators. For example, we have 

M*)M^) = ■ M%)Mx') ■ + Sf(x - x'), (28) 
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where the "abnormal" contraction S+ 3 (x — x') is specified below. The products ipp{x)%p a {x') 



and A^(x)A u (x') are expanded similarly, giving the abnormal contractions Si {x — x') and 
5^ u D + (x — x'), respectively. We have 

D + (x -x') = J dk B e- ik < x - x '\ (29) 

and 

iSf{x -x')= J dk F (}6±m F ) af3 e~ ik < x - x '\ (30) 
where the invariant phase space measure is 

Mb,f = -^6(k )6(k 2 - m 2 B F ) (31) 

with tub and mp are the boson and fermion masses. 
Now, using Dyson's parameterization (p2|), we find 

00 itn ( J\n [ m+n °° 



hFW = * E - } [ -^- II / *TiG(r t ) E^foryvgf*, (32 ) 

m % m\ n\ 1 iL , J I 
where the latter sum represents 



oo (F+F+B d 4 } 

e= e n (33) 

and A/^f s is the normal-ordered product of fields 

A^f 5 = : ^(Pi) ■ ■ ■ ^(Pf)HPf+i) ■ ■ ■ ^(Pf + f)A(Pf+f+i) • ■ • ^(Pf+f+b) : , (34) 

with 

iP(jp) = J d 4 xe ip - x ip{x) (35) 
= 2n6(p 2 - m 2 F ) £ {9(p°)u^\p)b a (p) + 0(-p°) v W(-p)dU-p)} , 

a 

■0 (p) = / d 4 xe~ ip ' x ip\x) and = / d A xe ipx A(x) = A'(—p). Evidently, we have 



oc 



E-M^AT M = J! / <*i r o / dte^'^^t-to) (36) 

M I *=° -oo J -oo 

xT { e ^ n+m tn +m ^( tn+m) . . . e -W" + i# /(tn+1 )} r {^o^) . . . e ^* 9jFf/ ( tn )| 

where A& = i(r'o+ri + ...+r/ c ). Assume that m > 1. Dyson showed that in the diagrammatic 
expansion based upon fl3"E| ) only connected graphs contribute, and therefore Atj%l must 
contain some L > 1 integrations over momenta ki corresponding to "abnormal" contractions 
between operators at times i > n and times tj, j < n. From the definitions above, we see 
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that these abnormal contractions will contribute Fourier components e~ lki '^ Xi ~ x ^ with the 
zeroth component kf > 0. We have then kf > mm(rriB, mp) for all L abnormal contractions. 

Next we make the shift of integration variables tj — > tj + t for j < n and consider Aij^l 
as a function of the variables Tj and tub^f- In calculating A4, we consider the Tj to be pure 
imaginary numbers in order to perform the time integrations, which have been extended to 
positive infinity. After these changes, the integration over t in (|36| ) yields a delta-function 

n+m n N 

M^ n (co, T) oc 6(co - ia £ T, - £ Wj - £ (37) 

i=ra+l j=0 Z=l 

where a>j is the algebraic sum of the zeroth components of the external momenta pk that 

F F+F+B m+n 

meet vertex j, so that u = Pk ~ S P° = J2 00 i- If we assume that mg and m F are 

fc=l fe=l i=0 

large enough, then 

L n+M n 

fcf > Lmin(mg, mj?) > uj — ia ^ Tj — J^cUj (38) 

Z=l j=n+l j=0 

so that ( |3~7|) cannot be satisfied and vanishes for m > 1. Thus for Re(Fj) = and 

itlb^f satisfying (^), M. may be found from (|36|) with m = 0, which then is simply an 
off-energy- shell Feynman amplitude. Dyson argued that considered as a complex function 
of the variables rriB,F and Tj, M. is analytic, and thus this representation of Ai is unique. 
So to obtain the intermediate Hamiltonian, one need simply calculate M. with the variables 
constrained as above (so that one can use generalized Feynman graphs), then analytically 
continue first mg^ from large values down to their physical values and then the Tj from 
pure imaginary back to pure real numbers (in practice, this analytic continuation is trivial). 
Finally, the integrals over the Tj in ([32] ) may be performed. 

C. Doubled Feynman graphs and spurions 

Dyson described a graphical method£] using what he called "doubled" Feynman dia- 
grams, for calculating intermediate representation operators from expressions such as (p6|). 
A doubled Feynman diagram is an ordinary connected n + 1-point Feynman graph in which 
a simply connected "tree" of (fermion and/or boson) propagator lines are drawn double. An 
example is given in Figure 1. The doubled lines account for the factors e aViti that occur at 
each point Xi and spoil energy conservation. If it is possible to draw more than one doubled 
Feynman graph from a given Feynman graph, one of these is chosen arbitrarily and the rest 
are ignored (they are not distinct, each one being obtained from any other by a suitable 
change in the internal momentum variables). Dyson then describes in detail how to shift 
the momenta of the doubled propagators consistently to account for the dependence of M. 
upon the Tj. The rules are somewhat involved. 



T Although, sadly, his papers have no figures. 
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We instead outline here an alternate but entirely equivalent (and hopefully more trans- 
parent) diagrammatic scheme that borrows from a formulation of perturbation theory due 
to Kadyshevsky.B! To do this, we express the step function in fl3"6|) as a complex integration 
through the familiar representation 

r du e iu ^~ U)) f dv e^ 1 ' 1 ^ 
J 2ni v — is J 2m v — is 

define the four- vector A M = crn^, where n 2 = 1 and n ■ x = x°; and write Hj(ti) = 
J d 3 XiHi(xi). Then ( |36D becomes 

E-^A/m = S m0 f ^-=?JL / die-*™**-"* {ft / d%) (40) 
{p} 16 J u=o J J 

xT {e (r °- l,/)X xo n I {x )e riX - Xi n I {x 1 ) ■ ■ ■ e r " Xxn H I {x n )} . 

From (|32|) and (0), we find the Fourier transform of the interaction part of the total 
intermediate representation momentum is 

dv —iTr 



K(p) = ^ h F(p) = A, £ L2L n / ^(ro / ^-^ (4i) 

n=0 U=0 ) Z/l ^ 4fc 

x J d 4 xe-^ p -^ )x > x [n / 

li=0 



X 



T {e (r °- w)A - a:o ^ 7 (xo)e riA - :Cl H / (x 1 ) • • • e r ™ A - x "ft 7 (x n )} . 



We expand the time-ordered product in terms of operator products fl34]) and perform the 
configuration space integrals in (IT), so that the resulting coefficient is a function of 

z/, the n + 1 variables Tj, the iV external momenta p^, and the total momentum difference 
p. These amplitudes will consist of integrations of m > n Feynman propagators depending 
on m internal momenta qj.- Shifting the space-time variables in (|4T| ) as Xi — > Xi + x, we see 
that .M is multiplied by an overall 5-function 

M?>(p,{r 4 },i/) = (2^\p-p°+p I )Mi P H{T},u). (42) 

This arises from integration over x. Here p° is the sum of outgoing and p 1 the sum of ingoing 
externa momenta. The remaining integrals over the Xi also give (^-functions, but the factors 
Tj and v spoil energy conservation at the vertices. For each vertex x-i we get a function 

(2vr) V 4 )(p? - p[ + q° - q\ - + 8 i0 u)X) (43) 

where pf is sum of the external momenta leaving Xj, q\ is the sum of the internal momenta 
entering Xj, and so forth. 

The amplitudes Ai^ are connected n + 1-point Feynman diagrams whose dependence 
upon v and the Aj may be accounted for by drawing an additional directed dashed line 
as in Figure 2, which is equivalent to Figure 1. This is very similar to the spurion line of 
Kadyshevsky's diagrammatic scheme.0 Although the rules for our scheme are quite different 
than those of Kadyshevsky's, we shall still call our dashed line a spurion line. We treat the 
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variables A M = — A'V and Af = iX^ J2)=o Tj as momenta and stipulate that the spurion line 
carries momentum as shown in Figure 2 — namely, the spurion line entering vertex i has 
momentum Aj, and that exiting vertex has momentum A. The spurion line touches each 
vertex exactly once, proceeding in decreasing order from vertex n to 0. A relabeling of the 
vertices 1 through n will not produce a distinct diagram since ([H]) is symmetric with respect 
to <-> Tj, i ^ j ^ 0, and so for a given choice of vertex we may choose any ordering 
of 1 through n and ignore all others, for they just sum to cancel the factor — y in fl4"T|). We 
may think of the spurion line as ordering the different vertices in time; however, it is only 
the order of vertex that need be distinguished from that of the others. 

With the spurion line considered as carrying momentum in the way described above, 
we have the usual four-momentum conservation at each vertex; and the evaluation of the 
amplitude M. proceeds as that of an ordinary Feynman graph. The spurion line serves 
only to shift the internal momenta so that the dependence on the Aj is now carried by the 
Feynman propagators — this is precisely the role of the doubled lines in Dyson's formulation. 
There will also always be an overall 5-function of the form 5^\p — A n + A), which expresses 
conservation of total three-momentum but not energy. The dependence of M. upon the 
variables r and v comes solely from this latter 5-function, and so the integrals over these 
variables in (f|l]) will yield the similarity function described earlier, a function that by choice 
of G(T ) falls off at least as rapidly as \ when the total energy difference u> = p° diverges. 

The most obvious advantage of Dyson's scheme is the apparent invariance of the expres- 
sion (|4T|). The Tj and v are dimensionless scalars, and covariance is broken only by A = an, 
which was chosen to have a special form. We note now that the choice n = (1, 0, 0, 0) above 
gives the equal-time formulation described here, but one could just as well choose n = n + 
with n 2 + = 0, n + ■ x = x + = x° + x 3 . This latter choice of A ensures that (41) yields the 
light-front formulation of the Hamiltonian, with x + the "light-front time." It is important, 
however, that divergences arising from the integrations in ( (4~i~l) are independent of A, for 
they occur when the internal momenta of our Feynman graphs diverge with respect to the 
variables Aj. Thus the determination of the counterterm structure of Hi is entirely equiv- 
alent to ordinary Feynman perturbation theoryjj] In particular, the usual rules for power 
counting, separation of divergences and multiplicative renormalizability apply, results that 
are by no means clear in other schemes for light-front variables. 



D. An example 

As an example we shall look at the quark-antiquark potential in light-front QCD. We 
do not discuss at length any of the many interesting, enticing and puzzling features of light- 
front QCD. For more details see Ref. [3] and the many references listed therein. Here we 
shall only outline the calculation and give the results. When formulated using light-front 



^ As with Feynman perturbation theory, the light-cone gauge choice A + = n + ■ A = is nonco- 
variant and necessitates the introduction of counterterms that depend on n± = (1,0,0, ±1). Such 
counterterms, of course, cannot depend on a. 
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variables! and light-cone gauge, half the degrees of freedom of QCD are constrained and 
one may develop an effective theory in which the gluons and quarks are two-component 
fields.il For some calculations it is more appropriate to use the original four- component 
constrained theory00 — for example, to demonstrate that the Slavnov- Taylor identities 
are fulfilled. 00 In general, however, calculations are much easier in the two-component 
theory, and we use just this theory here. The Feynman rules for the two-component effective 
theory are given in the appendix of Ref. [9]. We calculate the coefficient of the term 
A/"| 2 p|° = : ip* (ps)^ (— P2)i J (— pM^iPi) '■ m the intermediate Hamiltonian to second order in the 
undressed coupling g^. Perry" first showed using a discrete similarity transformation (with 
the similarity function / chosen to be a step function) that one obtains a logarithmically 
confining potential in the dressed Hamiltonian. ZhangQ later showed that a continuous 
similarity transformation such as those described here (specifically, the WG scheme of Ref. 
[3] with / a step function) yields a potential of the same form. Wilson and Robertsonlll 
gave general arguments (and caveats) for the appearance of such a confining potential in 
similarity transformed light-front QCD Hamiltonians. Here, we show how such a potential 
arises within Dyson's scheme. 

The graphs that contribute at second order are shown in Figure 3. Graph 3a comes 
from a term 7i\ (x) in the bare Hamiltonian that arises from eliminating the constrained 
components of the gluon field.§ This term is proportional to (g^) 2 , and in the spirit of our 
prior discussion, we take the term 7ij(x,x r ) appearing in the definition (|20D of Sd to be 
proportional to g 2 (X-(x — x')). The change to our momentum space Feynman rules caused by 
this modification is slight: we need only alter the way the spurion line carries the momenta 
T. If each vertex i in a given n+ 1-vertex diagram corresponds to a term in the Hamiltonian 
proportional to {go)™ 1 , then we have m = J27=o m i factors T^, 1 < j < rrii, corresponding 
to the functions g mi (X ■ (x — £«)). Then the diagrams are exactly the same, except for the 
redefinition Aj = iX J2)=o ^k=i ^j,k- The integration in (f41|) is now over 



The amplitudes corresponding to the graphs in Figure 3 are thus 

iM a = -Atg 2 B T a T a j-^-^(27r) 4 6^(p - A + A), (45) 



where p = p 3 + p 4 - p 2 - p\ and A = iX(T 01 + T 



02 ) 



^Light-front longitudinal momentum is p + = n + ■ p = p° + p 3 , light-front energy is p~ = n_ • p = 
p0 _ an( j transverse momentum is p\ = p\ i = 1, 2. The scalar product is p ■ x = ^p + x~ + 
^p~x + — p±_ • where the notation pj_ ■ x± means p l x l +p 2 x 2 . Occasionally, it convenient to 
define the light-front 3- vector p = (p + ,p±) and the product p ■ x = p + x~ — p± ■ x±. On-shell 
particles satisfy p~ = — L ^ — , from which it follows that transverse momenta scale as masses but 
longitudinal momenta scale differently. 
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^H — — m—*^ M) j-?nr (46) 




cr± ■ P2± — vm j , a± ■ Pa± + im 
a i — a 



2[p+] ^ x 2[p. 



(27r)V 4 )(j9-A 1 + A), 



where p is as above, A = i\T i, Ai = 2Ar n + A , and g = p 3 — pi + A — A x ; and, finally, 
A4 C is obtained from M. b by the substitutions A <-» — A 1; and A — > — A . In general, we 
need to regulate the singularities as the longitudinal momentum q + = n ■ q — > 0. As is 
usual for Feynman diagrams in light-front variables, we use the Mandelstam-Leibbrandt 
prescript iontJ 

(47) 



[q + ] q + q~ + is q + + isdt(q ) 
For the four-point interaction of ([45]), we interpret this prescription to be 



Pi -Pa f , f Pa -Pa ] 2 



[Pi" - P^] 2 I (Pi" - Ps")(Pi - Ps ) + *e j I (Pf - P3")(P4 - P2 



Z£ 



(48) 



The potential V A (pi,p 2 ,P3,P4) is proportional to the sum of these terms, Ai a + A4b + M. c 
integrated over the variables T^. We are interested in how this potential behaves at large 
separations, for which the terms proportional to [g + ] _1 will dominate the expressions in the 
brackets in Summing the three diagrams and relabelling the T, we have then for the 



long-range part of the potential 



iV LR 8®(p) = -4g 2 R T a T a 5^(p) f dT G(T ) f dT.G^) . (49) 

J J P — 1(J{Lq-\-Li) 



(q f q + q (q ) 2 q + q~ 



\ (<l +c l~ + i £ ) 2 1 +< 1~ — qj_ + ie (q + q~ + is) 2 q + q — q 2 _ + ie j ' 

where q~~ = p% — pi — iaTi and q~ = P2 ~ Pi + icr^i- For small q = p% — pi = P2 — Pa, 
the energy difference p~ can be ignored with respect to a. We calculate Vlr(x~ ,x±), the 
Fourier transform with respect to q of (^) above under this approximation. We need to be 
careful, however, for V^_r(0, 0) diverges. This is a result of well-known infrared singularities 
due to the choice of light-cone gauge, singularities that are known to cancel (in perturbation 
theory, at least) when evaluating gauge invariant quantities. 

Now, when using light-front variables, we need to take into account that transverse 
and longitudinal variables scale differently. Thus there is a transverse mass scale \i and a 
longitudinal momentum scale p. The energy scale is then a = fi 2 / p. It is natural to use 
dimensional regularization to handle both ultraviolet and infrared divergences in our scheme 
by taking the number of transverse dimensions to be d instead of 2. Then we must introduce 
an arbitrary (transverse) mass scale, which we take to be /1. This means the similarity flow 
is governed by the longitudinal scale p. That the couplings now run with two scales will 
have profound implications at higher orders, but this does not concern us at present. Setting 
p~ — > in (|49D, generalizing to d transverse dimensions and taking the Fourier transform, 
we find 
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OC 



„1 ..2 ~J y A-p yy 

Vlr(x-, Xi ) = f-rr / <ir G(r„) / dr 1 G(r 1 ) ?r -?-J r < 5 °) 



where 7 is the Euler-Mascheroni constant and E\(x) = —Ei(—x) is the exponential integral 
function. This gives a long-range logarithmic potential independent of the choice of G(T), 
which just affects its strength. Indeed, we can use the relation E\(z) — > —7 — In z for z small 
to see that 

Vlr ~ In p\x~ | for Ap\x~\ ^ Tin 2 x\] (51) 

otherwise, E\ is of order 1 and the In fi 2 x 2 ± term dominates. The divergent piece in Q5DD can 
be seen to be of infrared nature since it arises from a factor r(~ — 1), which is regulated by 
increasing d. 

We have dropped a short range term from the potential in going from (|9j) to (pUD, 
namely: 

gWT a T« 7 ^ 7^ 2r G(r )r 1 c?(r 1M . 



VW*.l) = ^ y ^r y dTi u v 17 < -icos- 1 ^ -P3±) (52) 



-/ cos 

Nj.||P3±l N±I|P4±I 

The dependence on the transverse momenta and in this term arises from the 
Mandelstam-Leibbrandt prescription for the infrared singularities. This potential is clearly 
small at large x±; however, at small x± it gives a — In p 2 x\ contribution that precisely can- 
cels the logarithmic term in (|50|) . Of course, to find the correct short distance form one has 
to also include the other terms from (pj) and (15) that were dropped in arriving at 



Finally, a simple choice of G that satisfies the above conditions is G(T) = 5(T — 1) + 
5'(V- 1). This yields 

V L r(x~,Xx) = I _ + 7 + 2* j + In ( ^ j ) , (53) 

where we have dropped yet another short range piece. Since the divergent constant in ( |53|} 
is infrared in nature, it does not get cancelled by a counterterm. This is important. Perryta 
has argued that this divergence is cancelled by a similar infrared divergence in the quark 
self-energy for color singlet states, but that for color nonsinglet states this cancellation does 
not occur. 



IV. TOWARD A NONPERTURBATIVE SOLUTION 

A. Dyson- Schwingerish Equations 

Within the WG and W schemes, there does not appear to be any great advantage in 
solving for the dressed Hamiltonian by using the perturbative expansion fllBl). Rather, we 
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can express the dressed Hamiltonian from the outset in terms of normal-ordered products 
of field operators as follows: 

HUt) = Y, h t{Pi,-,PF +F+ B)e ip0t M PFB = £^?-%}(*)> (54) 
to to 

where p = p° — p 1 , as before, and = h\ p J (2tt) 3 5^(p) — that is, total 3-momentum is 
conserved. Taking the Fourier transform of ([12]) — we only consider the W scheme here - 
we find 

oo 

H Ia (t) = J dt'f(t - + J da' J dt'dt"F au ,(t - t", t" - t')H Ia ,(t')H Ia ,(t"). (55) 

<j 

where f(t) is the Fourier transform of f{uj) and 

r±L^^m \ d ^m -WM) ,56) 

Now let us choose the similarity function motivated earlier, namely, f(x) = (1 + x 2 )^ 1 . Then 
the term in the brackets in (|56|) has a pole structure in u' that yields step functions in ti 
after integration. We find 



oo 



J2 ^2 



Hla (t) = °J dt'e-^Htf) -^J do'°—^- J dt'dfe-^'-^ (57) 

x ( e -|t-f| + e -^-*"l) [T{H Ia ,(t')H Ia ,(t")}-T{H Ia ,(f)H Ia ,(t")} • 

This is a very happy form. Indeed, when we substitute the expansion (5l[) in (0), we find 
that contractions of field operators resulting from normal-ordering the righthand side will 
only be of the Feynman type. Moreover, we need only consider the operator 



oc 



J2 ^2 



K a (t) = -j dt'e-^Hj(f) --J do'°-—^- J dt'df'e-^'-^ (58) 
[ e -«\t-t'\ + e-^'\)r{H Ial {t')H IfTl {t")} , 



since from the Hermiticity of Hi a we have 

H Ia (t) = K a (t) + Kt(t). (59) 

The calculation of K a involves only proper Feynman propagators, and one need only consider 
connected graphs since disconnected graphs will not contribute to Hj a (they are cancelled 
by identical terms from K\). 

We expand K a (t) as we did for H Ia {t) in (£§ — so that h\ p J = 4 p} + (k^Y, where 
kp = K( PP+F ,... ,Pi, — Pf+f+b, — Pf+f+i) — an d substitute in (|58|) . The amplitudes 
k^ so defined are to be solved by iteration. One can always organize such an iterative 
solution to reproduce the perturbative expansion of Section III, should one care to. Indeed, 
by expanding the dressed amplitudes in powers of the undressed amplitudes h\ p \ one 
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recovers the perturbative scheme of Section III. A. However, now there are other more in- 
teresting possibilities, for just as with Dyson-Schwinger equations one can use a variety of 
approximations to (|58| ) in order to sum up infinite classes of diagrams. We note that ( |58|) 



permits a Wick rotation to Euclidean space and thus the normal power counting applies in 
determining the counterterm structure. This does not seem to be true for all choices of the 
similarity function /. 

B. In the Dyson picture 

We can also develop Dyson-Schwinger-like equations for the intermediate Hamiltonian 
by introducing a new time scale T > into the integrals in (^3). We define 

S D (t,T) = T^e~ l ^-T dt ' Hl ^ t ' ) Y (60) 

so that 

to 



H?(t,T) = / dt Q g'{a(t-t )) J dt lQ {a{t - t x )) ■ ■ ■ (61) 

n=0 t _ T t _ T 

tn-l 

J dt n g(a(t - t n ))[#i(t n ), . . . [#j(*i),ffi(t )] • ■ ■] 

;-T 

t 

-a J dtog'iait-t^H^+i J dt l9 (a(t - ^[H^h), H?(t,t - h)], 



t-T 

t 



t-T t-T 



which is a linear integral equation. The usual intermediate Hamiltonian is then obtained 
from Hf(t,T) by taking the limit T — > oo. Hf(t,T) depends on t as an interaction rep- 
resentation operator — namely, Hf(t,T) = e' lHot Hf > (0,T)e~ lH ° t — so we expand Hf(t,T) 
in terms of normal-ordered products of field operators each at time t, with coefficents that 
depend on the positions of these operators and T: 

H?{t,T) = Y J h^{T)M {p] {t). (62) 
M 

N{p}{t) can be expressed as a normal-ordered product of field operators such as ip Pi (t) = 
ip(pi)e~ ip i t . Ko{t,T) is defined analogously to K a {t) above, and we find 

T r 

E^WrfW = / dt'J dTG{T)e-^' (63) 



p 



-i £ hf } h%\t')T{U {p} {t-t')U {pl} 
{p},{p'l 



The time-ordered product is expanded using Wick's Theorem in momentum space — for 
example, 
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T {^mlr(t')} =:1> q (t)i$tf): +i{2^\q - S F^e'^^. ( 64 ) 

We represent the product of all contractions — that is, the product of B q boson propagators 
D F (q i ), F g ipip^ propagators S^^Oto and F q ipty propagators — Sf( — <?fc)7o — as -M{g}> so 
that we have 

T{N {p] {t-t')N {pl] (t)} =5 NqNql £ \U( 2 ^ ) (4-m (65) 

xM {q} ■N {p / q} (t)N wm {t): e- ip0t '. 

The notation J\f{ p / q } means the N q = F q + F q + B q field operators with p^ = qi (those 
corresponding to the contractions) are dropped from the normal-ordered product A/| p }. Ex- 
pressing the amplitudes as 

kf{T) = dT 6 fcfr>(r) (66) 

J crl + ip" 

o 

and similarly for /i}^(T), after a bit of manipulation (|63|) becomes 

This equation has a simple graphical interpretation in analogy with the spurion diagrams 
of Section III.C. If we define the dressed amplitudes kffi(T } u) = 2tc5(p — iaT — av)kp (T) 
and h\ P p(r, v) = 2tc5(p° — iaT — av)lif^{T) and the undressed amplitude h\ p \T,u) = 
2n5(p° - iaT - av)G(T)h\ p} , then (|67]) becomes 

k^(T,u) = ^h\ p \T,u) (68) 

- .^}/tv^a> /ife /p, ^ } (r - ry - u') 



{,} F'F'B' v-v-ie 

If we draw a spurion line as in Figure 4, entering the undressed vertex with momentum 
irA M , carrying momentum {—v' + iT — iT^X^ from the undressed to the dressed vertex and 
exiting the last vertex with momentum — z/A M , then the total 4-momentum is conserved at 
each vertex. 

Figure 4 is a graphical representation of (|68|) for the A/^ amplitude. Undressed vertices 
are denoted by a dot, and dressed vertices (of type ko or ho) are denoted by a bubble. The 
resemblance to Dyson-Schwinger equations is evident. The dots indicate that we have not 
shown diagrams that modify the quark legs only (counterterm diagrams from the undressed 
Hamiltonian and a-dependent mass terms from the dressed Hamiltonian) as well as diagrams 
that involve dressed vertices with more than five total quark and gluon legs. The diagrams for 
fcjj are equivalent to those shown here, except the direction of the spurion line is reversed. 
Diagrammatically, it is easy to see how an iterative solution of these equations, starting 
with the undressed vertices as input, will reproduce the peturbative diagrams of Section 
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III.C. Of course, the diagrams are precisely equivalent only in the limit T — > oo. Note 
finally that because the amplitude k D (p 3 , — p 2 ; — p±,Pi] V) enters K D as the coefficient of 
A/" 220 (p3, — P2', — Pi, Pi), which is antisymmetric under either pi <-> — p 4 or p 3 <-> — p 2 and 
symmetric under (pi,—Pi) <-> G»3, — P2) plus complex conjugation, the diagrams in Fig. 4 
represent additional processes — for example, the last diagram shown represents processes 
where the undressed qg — > q vertex is attached to any of the four quark legs. 

The advantage of having an equation that is linear in the dressed vertices, as in the 
present case, is that the undressed Hamiltonian has a finite number of terms, whereas 
the dressed Hamiltonian can have terms with any number of legs. Thus a diagrammatic 
expansion for any given dressed vertex based upon (|58|), which is nonlinear in the dressed 
vertices, will contain an infinite number of graphs since there can be any number of internal 
lines connecting the two dressed vertices. In contrast, the diagrammatic expansion for any 
given dressed vertex based upon (^) contains a finite number of graphs, since the number 
of internal lines is limited by the undressed vertices. 



C. Another example 

As an example of this final formulation of the similarity transformation, let us set up 
the calculation of an approximation to the general quark-antiquark potential in light-front 
QCD. We start from the diagrams in Fig. 4. There will be similar diagrammatic equations 
for the dressed qg — ► q and qqg — > qq that appear in the equation for the dressed qq — > qq 
amplitude. We may approximate these latter equations as shown in Fig. 5. According 
to this approximation, we treat perturbatively vertices that change gluon number; in first 
approximation the dressed q —>■ qg is equivalent to the undressed vertex. Consistent with 
this approximation, we take the amplitude qqg — > qq to be as in Fig. 5. This corresponds 
to taking the input potential for the dressed qq — > qq amplitude to be as in Fig. 6. This 
yields the 0(gj l ) potential calculated in Section III.D. (For simplicity, we have neglected 
to draw the spurion lines in Figs. 5 and 6.) The equation for the dressed qq — > qq vertex 
in this approximation is shown in Fig. 7. This looks very much like the Bethe-Salpeter 
equation. Here again, we have ignored corrections to the external legs. Clearly, improving 
the approximations for the gluon- number- changing vertices corresponds to including higher 
order radiative corrections to the input potential of Fig. 6. 

While we shall not solve for this dressed amplitude here, it is perhaps instructive to set 
up the equations in some detail. From the above, the dressed vertices satisfy 

OO ■ rp 

h { rkt) = e** hm fdTf f-.^—h^(T, v) (69) 

T^oo J J ZlTl V — IE 


°° A 

= e^fdt'J dT J ge-^h%(T,u), 


with 

fc{g(r, v) = *r hf } (r, V )+ fdV E (70) 

{l} F'F'B' 
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X 



- ^—h { j P D ' q} (T\ v>) [M {q} ] 1 hf /p '^(T - T> - u>)} 



To ease our notational burden, let us call the dressed qq — > qq vertex = $; the dressed 
q — > qg vertex h} 1 ^ = 4>; and the dressed qqg — ► gg vertex /iff} — > \&. Undressed vertices will 
be denoted by the subscript — note that = 0. Now the approximations of Fig. 5 are 



~ rf>„(I»; (71) 



,<r, ^ijvj^j^. {^—ur, v>)s F{q w - r> - A (72) 

v ' 

+ — ^$(r', v')s F ( q )Mr - r> - v>) 

u' — IS 

The undressed quark-gluon vertex isi 

o (r, v) = G(T)M^) 4 t (4) (P3 -Pi-Qo- (*r + i/)A) (73) 
Mp^) = -9rT |2— ^ a ± -a ± ^ j. 

The undressed qq — ► gg vertex is 

$ (r, i/) = / dr G(r )dr 1 G(r 1 )5(r - r - ri)$ (27r) 4 <f( 4 V - (iT + ^)A) (74) 

$ (P3, -P2l -P4,Pl) = 4^ T ^ ] ^^ fF . 

We have seen these undressed vertices in the perturbative calculation of Sec. III.D. 
The next to the last graph in Fig. 4 is 

I dr 1% / (Cfi^— 4(r, WaW - r> - »>). (75) 

Using the approximation ([71]), this becomes 

$g(r, u) = -f dr G , (r )dr 1 G(r 1 )5(r - r - r x ) / ^ ^° . (76) 

J J Ziri v — v' — IE 

*4(-P4, ~P2, -q)D(-q)Mp 3 , Pl , q){2TxY5^\p - (iT + u)X), 

where q = P2 — Pa + (iTi + v')X. There is a similar contribution to $ from the term k ID . 
Indeed, it is identical after the replacement q — » q = p% — p% — (iT\ + v')X — this corresponds 
to reversing the direction of the spurion line. We write this as $q = $Q| 9 _>g. After integration 
over z/, the sum of these terms and the undressed vertex $0 gives the input potential of Fig. 
6, namely: 

v(t, v) = ar$ (r, v ) + $g(r, v) + $*(r, u) (77) 
= a J dr G(r ) J dr 1 c(r 1 )r 5(r - r - r 1 ) 

xF w (ir! + v)(2ir) 4 5^(p° - (iT + u)X), 
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where 

V {p} {u + iTi) = i» (g) + $ (g) - $>(-Pa, ~P2, -q)D{q)Mp 3 , Pl , q) (78) 

-0o(-P4, ~P2, -q)D(q)(f)o(p 3 ,p u q), 



with g = p 2 — Pi + [y + ^r^A and q = p 3 — p\ — [y + ir^A. Substituting this into (|69]), we 
recover the order potential of Sec. III.D, as indeed we must. Thus the first and fourth 
terms on the righthand side of the equation in Fig. 4 (plus their conjugate graphs) sum to 
give the first term on the righthand side in Fig. 7. 
The last graph in Fig. 4 is 

I dr ' lik l$?!T=h&®'< ^to*(- 9 ')*(r -!>-✓)■ (*» 

Using the approximation (0), the integrand in ([79]) becomes 



dv" r d^a" 1 
' 2vri J (2vr) 4 i/ - i/ - p" - is y 1 



x 



, ^-^-«(r', z/) 1 D{q)S F {-q')^^\Y\ u")S F (q")<S>(r - V - T", v - u' - u") 



plus a similar term in which the gluon couples leg — p 2 to p 3 . Other terms are either ignored 
because they only correct the external legs or else not distinct because of the symmetries 
of the amplitude. The term given in (|8~0D may be summed with the second graph on the 
righthand side of the equation in Fig. 4 to give 

[dv' r d\ r d A q' , , , S F (q)S F (-q') 

dT I^JJ^yJ j^ v ^-^ {T > u) u-i/-ie lw *• " i -* ,T " r ' " ~~ v ] (81) 

This is the second graph on the righthand side in Fig. 7. The remaining graphs in Fig. 7 
result from (|H]) by complex conjugation (reversing the spurion line) and by the substitution 
—P2 — Pi- We have, finally, 



^P3,-P2;-p4,pi 



+ ®q-p2 ~q',Pl\ l > Z> ) , _ ■ V P3-q'--P4,q{ L ~ 1 J V 

Is oc. 



-Vq,q'; P 4,pA L > v ) ~ _ ^ _ ^. ^PB^g'.gl 1 _i > ^ ~ ^ J 

- 4 *?,3';P4,Pil i > ^ J — ^ _ i£ — v Ps,P2\q',q\\ ~ 1 -> v ~ v ))■ 

It would be interesting to approximate V further by just taking the long range (small q) 
part calculated in Sec. III.D and substituting this in (|82|). The divergence in V will be 
cancelled by divergences in the self-energy corrections to the external legs (see Sec. III.D) 
- corrections which we have ignored here but which must be included if one wants to solve 
this Bethe-Salpeter-like equation in a consistent manner. This paper is long enough already, 
and therefore such a calculation will not be undertaken here. 
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V. CLOSING ARGUMENTS 



In this paper we have developed a similarity renormalization framework for the study 
of strongly interacting field theories. The basic view of this approach is that Hamiltonian 
methods will prove most fruitful for the solution of such theories. The generic theory for 
which our framework has been developed has interactions that are weak asymptotically but 
grow large enough at low energy to invalidate perturbation theory. QCD is of course the 
theory that we particularly have in mind. The basic idea behind the similarity renormal- 
ization scheme is that the Hamiltonian of such a theory can be transformed into that of 
an effective many-body theory where states that differ greatly in free energy are essentially 
uncoupled. This resulting many-body Hamiltonian can then be solved with recourse to the 
many available approximate methods such as variational approaches, trial wave functions, 
iterative techniques and numerical basis function approaches. The similarity transformation 
achieves this goal by smoothly eliminating interactions between states whose free energies 
differ by an amount greater than an arbitrary scale a introduced by the transformation. The 
result of this elimination is to create new interactions between states that differ in energy by 
an amount less than a. The assumption is that from these new interactions one can extract 
potentials that describe the relevant physics. 

Clearly, no physical result can depend on the renormalization scale a. However, the 
calculations can only be approximately carried out for most systems of interest, and so one 
will need to choose a value of a that makes the approximations as valid as possible. Com- 
putations in this framework involve two steps: first, the similarity transformed Hamiltonian 
is calculated; second, this Hamiltonian is solved. The first step is most easily done per- 
turbatively, which requires that a be large since only at high energy is the coupling small 
(by assumption). The second step is most easily done when as few states as possible are 
coupled strongly, which requires that a be as small as possible. Inevitably, one must make 
a compromise and choose a in some intermediate range. If the Dyson-Schwinger-like equa- 
tions developed in the last section can be implemented, it may be possible to go beyond the 
perturbative calculation of the dressed Hamiltonian so that a can be lowered to values for 
which perturbation theory begins to break down. However, in most cases it will probably 
be better to keep a large enough that the Hamiltonian can be dressed peturbatively, since 
the tools available for the second step are generally more powerful. 

That the elimination of couplings between states separated by large energy differences 
turns a field theory into a many-body theory is clear if the free particles are massive. If 
a is lowered to just above the particle masses, then creating or destroying a particle will 
cause an energy difference of order a. This means interactions that change particle number 
will be suppressed. If the particles are light or massless, however, as in QCD, it does 
not follow that the similarity transformation described here will suppress interactions that 
change particle number. Indeed, one might consider it better to have the similarity function 
depend on differences in particle number rather than on differences in energy. However, 
such a transformation is unlikely to be feasible in perturbation theory, for one no longer has 
a means of eliminating only high energy degrees of freedom. Moreover, one would lose the 
advantage that small energy denominators are avoided. Instead, one can use the similarity 
transformation in a light-front formulation of a theory such as QCD. We shall argue in the 
following that this change of formulation of the theory will provide a means of turning even 
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QCD into a many-body problem. 

The use of light-front coordinates in field theory is not new (see Ref. [3] for a review and 
an exhausting — if perhaps not quite exhaustive — list of references). Such calculations 
are considered difficult, partly because they are unfamiliar, but also because they really 
are difficult. One of the major difficulties with the use of light-front coordinates is that 
Lorentz covariance is no longer manifest once one chooses a particular direction in defining 
the light-front time. Another problem is that calculations are only feasible for a particular 
gauge choice, the light-cone gauge. These problems severely complicate the renormalization 
of divergences, and this has slowed advance in light-front computations. The fact that in 
our similarity renormalization scheme the calculation of the dressed Hamiltonian can be 
expressed in terms of Feynman diagrams is therefore significant. It means that divergences 
can be handled with dimensional regularization, the infrared singularities can be handled 
with the Mandelstam-Leibbrandt prescription, and the counterterm structure is then that 
of covariant formulations. In particular, renormalization of divergences for QCD in equal- 
time coordinates and the light-cone gauge is well understood (see, for example, Refs. [10- 
11]), and we can rely on such prior work in determining the dressed Hamiltonian in the 
intermediate representation. There are of course other, more technical difficulties resulting 
from the choice of light-front coordinates (for example, rotations about the transverse axes 
are no longer kinematical); but these should not be prohibitive. Indeed, there are also some 
technical advantages to the use of light-front as compared to equal-time coordinates (for 
example, now boosts are kinematical). 

The essential reasons for using light-front coordinates follow from the free particle dis- 
persion relation 

- p 2 ± + m 2 

P = r • 83 

p+ 

This implies that all on-shell particles (and antiparticles) have have longitudinal momentum 
p + > and all physical states have total longitudinal momentum P + > 0. The vacuum 
has P + = and is therefore built only from particles with p + = 0. From ( |83| ) we see that 
such particles have infinite (light-front) energy unless p\ = and m 2 = 0. The similarity 
transformation will decouple the infinite energy states from the low-energy physics, and the 
computation of the vacuum state is then greatly simplified. Indeed, if we give gluons a small 
mass, all particles with p + = have infinite energy, and the vacuum state of the dressed 
Hamiltonian is just the trivial free vacuum state. Likewise, in building any bound state 
of the dressed Hamiltonian, we can rely on the fact that the total longitudinal momentum 
P + must be built up from constituent longitudinal momenta pf > 0. If the energy of each 
constituent is to be less than a, we must have pf > {p]_ + m 2 )/a for all pf, and thus in 
effect we have a bound on the number of particles that can make up a state. Again, the 
only particles that escape this lower bound on the longitudinal momentum are (massless) 
gluons with p\ ~ 0, but one can hope that such highly transversely localized particles will 
have little contribution to low-energy states. In any case, it is clear that effects associated 
with the quark-gluon sea in the equal-time formulation of QCD will be largely, if not totally, 
replaced by interactions in the dressed Hamiltonian after a similarity transformation in the 
light-front formulation. 

But the most interesting result of the light-front formulation is that longitudinal and 
transverse variables scale differently (this has been repeatedly emphasized by Wilsonlll). 
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This can be seen from the dispersion relation ( |S3|) and has been discussed in Sec. III.D. 
Thus the similarity renormalization of light-front Hamiltonians depends on the energy scale 
a through both a mass scale p and a longitudinal momentum scale p. The mass scale p 
may be taken to be that introduced by dimensional regularization. The QCD coupling 
runs with this scale, and p must be taken large to keep Qr small enough for perturbation 
theory to remain valid. However, the similarity function / runs with the ratio a = p 2 /p; 
and thus, even though the mass scale p be large, we can suppress couplings between states 
with energy differences above a given value by making the longitudinal scale p large enough. 
We have seen that light-front energy increases with particle number; and therefore at large 
p the similarity transformation will suppress number changing interactions in the dressed 
Hamiltonian, thus turning the solution of QCD bound states into a many-body problem. 
That the mass scale p may be nevertheless kept large provides the possibility that the trans- 
formation may be accurately calculated in perturbation theory. That Dyson's intermediate 
representation is a particular similarity renormalization scheme that may be expressed in 
terms of ordinary Feynman graphs means that perturbative light-front QCD computations 
are feasible. These computations should yield a dressed Hamiltonian that is considerably 
closer to strong-interaction phenomenology than the usual free QCD Hamiltonian. That we 
find a confining potential already at tree level may be taken as a positive indication of the 
utility of this approach. 
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FIGURES 



FIG. 1. Doubled Feynman graph 




FIG. 2. Spurion graph corresponding to Figure |l] 




FIG. 4. Diagrammatic equation for general qq — > qq vertex 




FIG. 5. Approximations for q — > qg and qqg — > (7(7 vertices 




FIG. 6. Definition of input potential for approximation to — > qq vertex 




